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We derive all relevant equations needed for constructing a global general relativistic 
model of advectively cooled, very hot, optically thin accretion disks around black holes 
and present solutions which describe advection dominated flows in the gravitational 
! field of a Kerr black hole. 
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1. ADVECTION DOMINATED BLACK HOLE ACCRETION DISKS 

6 

in ' Recently, solutions corresponding to a new type of hot, optically thin accretion disks cooled 

^ ■ mostly by advection has been discovered (see Abramowicz &: Lasota 1996, Abramowicz 1996, 

and Narayan 1996 for reviews). The luminosity, L, of these disks is much lower than that of the 
standard disks with the same rate of mass supply M. The radiative cooling efficiency, defined 
as r] = L/McP, is typically less than 0.01 in the newly discovered disks, while in the case of the 
I standard disks it is typically ry ^ 0.1. Thus, the advection dominated accretion flows (ADAFs) 

hide the evidence that they are consuming a lot of matter — for this reason Abramowicz & 
Lasota (1996) called them "secret guzzlers". A typical secret guzzler appears to be very hot and 
underluminous. 
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Lasota, Narayan and collaborators have discussed properties of the secret guzzler models and 
compared them with observed properties of sources which arc believed to contain accreting black 
holes and appear both undcrluminous and very hot, like Sgr A* (Narayan et al. 1995), soft X-ray 
transients (Narayan et al. 1996b, Lasota et al. 1996a), NGC 4258 (Lasota et al. 1996b). These 
detailed discussions, based on calculations of electromagnetic spectra, convincingly demonstrated 
that secret guzzlers may explain these particular types of sources, with no need for ad hoc 
hypothesis about additional optically thin components in the model, or some unusual properties 
of the viscosity. 

Models of secret guzzlers are described by similar equations as the well-known slim accretion 
disks (Abramowicz et al. 1988). The only difference is in equations for the radiative processes: 
secret guzzlers are hot and optically thin and the ion temperature is, most probably, much higher 
than that of the electrons (see e.g. Narayan 1996) and slim disks are cool and optically thick. 
All the previously calculated models of slim-disks and secret guzzlers are not fully satisfactory 
at least in one respect: they assume an approximate model for the gravity of the black hole. It 
is either just the plain Newton's potential of a point mass M, $ = —GM/r, or the Paczyhski's 
pseudo-potential $ = -GM/{r - re), = 2GM/c^ (Paczyriski & Wiita, 1980, hereafter PW). 
No general relativistic solution has ever been obtained. 

The PW paradigm is an excellent approximation for the calculations of the hydrodynamical 
properties of accretion disk around non-rotating black holes. It is both very simple and very 
accurate. However, it has two shortcomings: (1) it does not describe the relativistic effects in 
the light propagation, and (2) it does not include the effect of the inertial frame dragging, an 
important property of rotating black holes. For this reason, at the TAD2 meeting at Garching 
(see Duschl at al. 1994) one of us suggested that the time has come to construct solutions for slim 
disks that would be fully consistent with general relativity. The first step should be a derivation of 
the relevant equations. Ramesh Narayan argued that it should be the job of the proposer so this 
was done in Lasota (1994; thereafter L94). 

Wc repeat this calculations here, correcting some computational errors of L94 and using a 
slightly more practical notation. The general assumption and definition are presented in Section 2. 
In Section 3 wc give the set of equations describing a general relativistic model of a hot accretion 
flow. The detailed microphysics that describes the radiative cooling is briefed in Section 4 and 
5. After a discussion of the transonic character of the flow in Section 6, numerical solutions of 
advection dominated flows in the gravitational field of a Kerr black hole are presented in Section 
7. Section 8 contains the conclusions. 
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2. ASSUMPTIONS AND DEFINITIONS 



We use geometrical units here that are hnked to the physical units for length, time and mass 

by, 



length \ / length 
in physical units J \m. geometrical units 



time \ 1 / time 
in physical units / c V in geometrical units / ' 



mass \ f mass 



in physical units / G v in geometrical units 



(1) 
(2) 
(3) 



We use the Boyer-Lindquist spherical coordinates i, r, 9, ip to describe the Kerr black hole 
metric. The cylindrical vertical coordinate z = cos is defined very close to the equatorial plane, 
2 = 0. The metric of the Kerr black hole on the equatorial plane, accurate up to the {z/r)^ terms, 
takes the form given by Novikov & Thorne (1973, hereafter NT), 



ds^ = —df + ^(dip- Lodtf + — dr^ + dz"^ (4) 



2Mar 

A = r2-2Mr + a^ A = + r'^a^ + 2Mra^ , oj = ^ . (5) 

Here M is the mass, and a the total specific angular of the Kerr black hole. The metric (Q) 
was also used in L94. Note, that the determinant of the metric tensor corresponding to @) is 
9 = -r"^- 

There are two Killing vectors in the geometry (^), connected to the time symmetry (metric is 
independent of t), axial symmetry (metric is independent on ip) 

ri^ = 5\,y e = 5\^y (6) 

where (j.-^ is the Kronecker delta. Using Killing vectors (|6|) one defines some useful scalar 
functions: angular velocity of the dragging of inertial frames a;, gravitational potential <I>, and 
radius of gyration R, 

^ = -||, e-2* = (w)-u;2(eO, R' = me'''- (7) 

In the Boyer-Lindquist coordinates the scalar products of the Killing vectors (P) are given by 
the components of the metric, 

(w) = 9tu = 9t^, {^0 = 9w^ (8) 
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and therefore quantities defined by (|^ can by explicitly written down in terms of the 
Boyer-Lindquist coordinates: 

r^A' A ■ ^ ' 

The black hole surface (event horizon) is at 

r+ = M + (^M^ - a'^y^'^ . (10) 

The unit timelike vector 



{v' + ^e), (11) 



is orthogonal to the space- like surfaces t =const. It corresponds to the four- velocity of the local 
inertial observer or ZAMO, i.e. Zero Angular Momentum Observers (Bardeen 1973). 

The four velocity of matter has components u*, u'^, vJ' , 

u^ = u'5\,^+un\^^ + u^5\,y (12) 

One defines the angular velocity with respect to the stationary observer, and the angular 
velocity with respect to the local inertial observer by, 

J7 = ^, Q. = 9.-uj, (13) 
It* 

The angular frequencies of the corotating (+) and counterrotating (-) Keplerian orbits are 

, Ml/2 

= -7^, (14) 

and the Keplerian specific angular momentum is given by 



Ml/2 ^^2 ^ 2aM^l^r^l^ + 



= ± ^ T79- (15) 



r 



3/4 (^3/2 _ 3MrV2 ± 2aMi/ 



The Keplerian angular momentum has a minimum at the marginally stable orbit 



r-, = M{3 + Z2T[(3-Zi)(3 + Zi + 2Z2)]i/2}, 



1/3 



Zi = l + a'/NP) '' ' [(1 + a/Mfl^ + (1 - a/Mf'^ 

1 /2 

Z2 = (Sa'^/M^ + Zf'^ . (16) 
In the reference frame of the local inertial observer the four velocity takes the form, 

n^=^(n^ + >)r^(^)+^Mry . (17) 
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The vectors t^^^^ and r*^^,^ are the unit vectors in the coordinate directions (p and r. The Lorentz 
gamma factor 7 equals, 

7= ; ^ (18) 

^1 _ (yMf _ („(r))2 

The relation between the Boyer-Lindquist and the physical velocity component in the azimuthal 
direction is, 

v^'^^ = M. (19) 
It will be convenient to use the (rescaled) radial velocity component V defined by the formula, 

-^—=^yir)=urgl/^. (20) 

The Lorentz gamma factor may be written as, 



7^ 



i-mO (i-fO' ^^^^ 



and from this follows a simple expression for V in terms of the velocity components measured in 
the frame of the local inertial observer, 

V = , = , , , ■ (22) 

Thus, V is the radial velocity of the fluid as measured by an observer at fixed r who corotates 
with the fluid. In the notation used in L94 (cf. his equation [17]), 



(r) 



(23) 

L94 ^ ' 



Although a different quantity could have been chosen as the definition of the "radial velocity" , 
only V has the three very convenient properties, all guaranteed by its definition: (i) everywhere in 
the flow \V\ < 1, (ii) on the horizon \ V\ = 1, (iii) at the sonic point \V\ ~ Cg, with being the 
local sound speed. To see that the first property holds, let us define = u''u,f = u''u'' g„- > 0. 
Then, one has V'^ = V"^ /{I + V"^) < 1. Writing V = 1/ + A/{r'^u^u^') demonstrates property 
(ii). Property (iii) of V will be proved in Section 6. Other possible choices of the "radial velocity" 
are not that convenient. For example, the radial velocity u = lu^], chosen by Shapiro & Teukolsky, 
(1983) has none of these properties and would be less convenient in the present context. 

The stress-energy tensor T*'^ of the matter in the disk is given by, 

T^*^ = {e+p) u'v}" +pg''' + S'^ + u^q' + u'q'', (24) 
where e is the total energy density, p is the pressure, 



Sik = i^pcrik, 



(25) 
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is the viscous stress tensor, p is the ress mass density and is the radiative energy flux. In the 
last equation v is the kinematic viscosity coefficient and 0",^ is the shear tensor of the velocity 
field. From the first law of thermodynamic it follows that 

de = ^^dp + pTdS, (26) 

where T is the temperature and S is the entropy per unit mass. Note, that in the physical units 
that e = pc^ + n, where 11 is the internal energy. For non-relativistic fluids, 11 <C pc^ and p <^ pc^, 
and therefore 

e + p « c'^p. (27) 

We shall use this approximation (in geometrical units e + p ~ p) in all our calculations. This 
approximation does not automatically ensure that the sound speed is below c, and one should 
check this a posteriori when models are constructed. We write the first law of thermodynamics in 
the form: 

dU = -p (-] + TdS (28) 



where U = H/p. 



3. SLIM DISK EQUATIONS IN KERR GEOMETRY 

It is convenient to write the final form of all the slim disk equations at the equatorial plane, 
z = 0. Only these equations which do not refer to the vertical structure could be derived directly 
from the quantities at the equatorial plane with no further approximations. All other equations 
are approximated — either by expansion in terms of the relative disk thickness H/r, or by vertical 
averaging. The thickness of the disk H is defined as in equation ( |4^ ) . 



3.1. Mass conservation equation 



From general equation of mass conservation, 

V' (pui) = 0, (29) 

and definition of the surface density S, 

r+Hir) 

S = / p(r, z)dz PS 2Hp, (30) 

J-H(r) 



we derive the mass conservation equation, 

M = -27rA^/^S , ^ (31) 
It is identical to that derived by NT. In L94 the numerical factor on the right hand side is — 47r. 
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3.2. Equation of angular momentum conservation 

From the general form of the angular momentum conservation, 

Vfc (r^'^i) = 0, (32) 

we derive, after some algebra, 

M dC I d .^i^d^/'^-f^ dVt 



2'Kr dr r dr \ dr 
where F~ = is the vertical flux of radiation, and 



+ rx: \-f~c = 0, (33) 



/ ^3/2 \ ^ 

C^-{uO = -u, = ^\-^\n, (34) 



is the specific (per unit mass) angular momentum. Equation ( |33| ) differs from that derived in L94. 
In addition to a trivial misprint in the L94 equation, there is also a more significant difference 
— we keep the term F~ C that was rejected from the final form of the L94 equation. This term 
represents angular momentum losses through radiation. Although it was always fully recognized 
that angular momentum may be lost this way, it has been argued that this term must be very 
small. Rejection of this term would enormously simplify numerical calculations, because it is 
obvious that with C = equation (|3^) can be trivially integrated. 



-(£-£„) = -E..4'/^^-, (35) 

where Cq is the specific angular momentum of matter at the horizon (A = 0). In the numerical 
scheme for integrating the slim Kerr equations (with F~ C assumed to be zero) the quantity Cq 
plays an important role: it is the eigenvalue of the solutions that passes regularly through the 
sonic point. 

Several authors have pointed out that in some astrophysically relevant situations the F~ C 
cannot be put to zero, because it could be a significant part (up to 0.4) of the angular momentum 
balance (see Lamb, 1996 for discussion and most recent references). 



3.3. Equation of momentum conservation 



From the r-component of the equation VjT*^ = one derives 

V dV A IdP 



1 — V"^ dr r S dr ' 



(36) 
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where P = 2Hp is the vertically integrated pressure and 



MA {n-nt)in-n,) 



L94 version of the momentum conservations equation contains few trivial misprints. 



Using equation (^) one shows that (|36| ) has the correct limit corresponding to the spherical 
accretion (il = 0) on the non-rotating (a = 0) black hole: in this limit it becomes identical with 
equation (G.7) on page 569 of Shapiro & Teukolsky (1983). 



3.4. Equation of energy conservation 

From the general form of the energy conservation 

V. (rSfe) = 0, (38) 

and the first law of thermodynamics, 

the energy equation can be written in general as 

^adv^^+_^-^ (40) 

where 

F--^4j.'f?)^ (41) 



'•^ \ dr 

is the surface heat generation rate, F~ is the radiative cooling flux (both surface) which is 
discussed in Section 5, and F'^'^'^ is the advective cooling rate due to the radial motion of the gas. 
It is expressed as 

^adv ^ T— = T— . (42) 

y/l - y2 r dr 27rr dr 



3.5. Equation of vertical balance of forces 



One should calculate the condition for the hydrostatic vertical equilibrium in the frame that 
comoves with matter. NT have approximated this frame by the one which moves along circular 
Keplerian orbits at the equatorial plane. It is known that in the innermost part of the disk the 
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motion of matter differs considerably from the circular Keplerian motion. For this reason, one 
needs to repeat the NT calculations in the comoving frame of matter. [] 

Like NT, we demand that the vertical pressure force should be equal to the vertical component 
of the tidal gravitational force, 

-p{^)^ = HR',,iTl^g^,T%^u^u^ « H^'R^^,,rl)g^,T%^n^n^ = ^T^i?^ (])(.)(,) , (43) 

however we demand this to happen in the comoving frame of the fluid (note the fluid's four 
velocity u^u'' multyplying the Riemann tensor R^jj^i), while NT demand that equality of forces 
should occur in the frame rotating with the Keplerian velocity (NT have w'p^v}^ multiplying the 
Riemann tensor). In the formula r'^^^ is the unit vector in the z-direction. From this we derive the 
final formula, 

P _ 2^ 

that should be compared with NT formula (5.7.2). The difference is that in our formula the 
redshift factor 7 is given by (^), which corresponds to the motion of matter and not to a fictitious 
fluid with strictly Keplerian rotation. The NT formula is singular at the location of the photon 
circular trajectory, r = rph, and our formula at the horizon. L94 gives the same formula for the 
vertical equilibrium. 

The singularity of Eq. (^) on the horizon is an artifact of approximations used in Eq. (p3[). 
We assumed that the flow is vertical hydrostatic equilibrium, in the ZAMO frame, down to the 
black hole surface. This requirement implies an infinite boost at the horizon. In reality the flow 
near the horizon is supersonic and therefore practically in a free fall. Eq. (^) should be completed 
by a term (dv^ /dr) which, in the supersonic flow, will dominate over the pressure gradient 
term. As a result Eq. (|4^ ) will become the (non-singular) geodesic deviation equation. However, 
because we are adopting here a ID approach in which 'vertical' quantities are averaged it is not 
possible to include Vz terms in the equations. Since the 'singularity' appears only in the supersonic 
part of the transonic flow and since the reason for its appearance is clear, it has no influence on 
the physically important properties of the accretion flow (see also discussion in Section 6). 



3.6. Viscosity prescription 



^This reason diflers from that discussed recently by RifFert & Herold (1995, Ap.J., 450, 508) who have also re- 
derived some of the NT formulae. They have argued that for some of these derivations, the Kerr metric must be 
accurate to second order in z/r. Although they have kept the (z/r)^ terms in geometry, they have rejected physical 
quantities of this order: radial pressure gradient, radial advective cooling, etc., which is not a self-consistent scheme. 



(H + a^Y + 2Aa^ 



(r^ -|- a 



2\2 



Aa2 



(44) 
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The standard assumption for the viscosity coefficient is, 

V = {2/3)aCsH, (45) 
where Cg = \ppjp is the isothermal sound speed. 

4. THERMODYNAMICAL RELATIONS 

The equation of state can be expressed in the form: 

p = PrA pTjH pTe + —, (46) 

Hi /ig 247r 

where pr is the radiation pressure, TZ is the gas constant, and fie. the mean molecular weights of 
ions and electrons respectively, Tj, and ion and electron temperatures, a the radiation constant, 
and B the intensity of a isotropically tangled magnetic field, includes the radiation, gas and 
magnetic pressures. The radiation pressure pr, the gas pressure pg, and the magnetic pressure 
Pm correspond respectively to the first term, the second and third terms, and the last term in 



equation (46). 



The mean molecular weights of ions and electrons can be well approximated by: 

where X is the relative mass abundance of hydrogen and Y that of helium. We may define a 
temperature as 

T = J^ + '^), (48) 

where 

_ / 1 1 \ 2 
^-[j^^yj ^ 1 + 3X + l/2y ^^^^ 

is the mean molecular weight in the standard approximation (Cox & Giuli 1968). In the case of a 
one-temperature gas (Tj = Te), one has T = Ti = T^- For an optically thick gas, pr = ^aT^. 

Until the end of this section we will assume that p^ = ^aT^. As noted by Narayan &; Yi 
(1995b; hereafter NY), the frozen-in magnetic field pressure pm ~ ~ p'^^^, therefore we may 
write the internal energy as 

U = ^+ ^, +e.pV3^ (50) 

p pruuijg - 1) 
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where Cq is a constant {pm = l/Sco/o^''^) and 7^ is the ratio of the specific heats of the gas. We 
define 

(51) 



/3 = ^, /3™ 



P Pg + Pm ' ' 3/3, 

From equations (^6|) and ( ^0|) one obtains the fohowing formulae (see e.g. Cox &: Giuh 1968) for 
the specific heat at constant volume: 



7^ 



12(l-/?//?^)(7g-l)+/? ' 
/3 



4-3/3* p 



and the adiabatic indices: 



r. -1 



(4-3/3*)(7g-l) 
12(1 -/?//?„) (7^-1)+/? 



(52) 
(53) 



Ti =/3* + (4-3/3*)(r3-l) 



(54) 



The ratio of specific heats is 7 = Cp/cy = Ti//?. For j3 = j3m we have Fa = 7^ and 
Fi = (4 — /3)/3 + /3(7g — 1). For an equipartition magnetic field {(3 = 0.5) one gets Fi = 1.5 and 
for (3 = 0.95, Fi = 1.65 (here we have used 7^ = 5/3). Our formula for Fi is different from that 
used by NY but the numerical values differ by less than 5%. One expects (3m ~ 0.5 — 1. Since 



dr 



cy 



dlnT 



dr 



(ra - 1) 



dlnS d\nH 



dr 



dr 



(55) 



In the numerical scheme based on Chen & Taam (1993) (see also Abramowicz et al 1995) the 
advective flux is written in the form: 



F 



adv 



M p 



27rr^ p 



where 



3(3*dlnT ^^,,dlnS 
+ (4 -3/3*' 



(56) 



(57) 



F3 — 1 d In r d In r 

Here the term oc din H/ din r is neglected. Since no rigorous vertical averaging procedure exists, 
the presence or not of the din H/dlnr - type terms in this (and other) equation may be decided 
only by comparison with 2D calculations (see e.g. Narayan & Yi 1995a and also Chen et al 1995 
for a discussion of this point). 

For /3m = 1 we recover the formulae used in Abramowicz et al (1995). The formulae derived 
in this section are valid for the optically thin case r = if one assumes /3 = Pm ■ 



5. RADIATIVE COOLING 
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We have to complete the set of equations by specifying the physical processes which will be 
involved in the radiative cooling. In our opinion the most convenient and general description of 
cooling processes have been presented by NY. They consider a two-temperature plasma cooled 
by synchrotron radiation, inverse Compton process and bremsstrahlung emission. They neglect 
electron-positron pair creation and annihilation but as shown by Bjornsson et al (1996) and 
Kusunose &; Mineshige (1996) this is justified in most cases of interest. Below, for the reader's 
convenience we will recall the formulae used by NY. 



5.1. Heating of electrons by ions 



NY use the formula of Stepney & Guilbert (1983) which has been modified in order to account 
for a X = 0.75 and Y = 0.25 composition (the effective molecular weight of ions /Xj = 1.23, of 
electrons = 1.14). Here we give the general formula in which the factor rii in the Stepney &; 
Guilbert (1983) formula is replaced by, 



n 



3 "'3' 



(58) 



where Zj and rij are respectively the charge and the number density of jth species. The Coulomb 
collisions transfer energy from (hotter) ions to electrons at a volume transfer rate 



r = 



3 TUp 



2mi"^"^^''K2(l/ee)i^2(lM) 



In A 



+ 1 



+ 2Ko 



erg cm ^ s 



(59) 



where the K^s are modified Bessel functions InA 20 is the Coulomb logarithm and the 
dimensionless electron and ion temperatures are defined by 

9e = A;Te/meC^ 9i = kTi/rUeC^. 



(60) 



5.2. Bremsstrahlung Cooling 



The bremsstrahlung includes emission from both ion-electron and electron-electron collisions 
(Stepney & Guilbert 1983; Svensson 1982; NY) 

/br = /ei+/ee- (61) 

The ion-electron bremsstrahlung cooling is given by 

f~=neriaTcafme(?Ff,i(6e) (62) 
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where a/ = 1/137 is the fine structure constant and 



F,i{e,) = ^[-^j (1 + 1.7810,^-^^), 0e<l, 
90 

= — pn(1.1236le + 0.48) + 1.5], Oe > 1. (63) 
zvr 

In the original formula quoted by Stepney & Guilbert (1983) there is a number 0.42 instead of 
0.48 (see NY). 



In the non-relativistic approximation one gets from Eq. (62) 

/- = 1.57 X lO^^'^ngnTg^/^ erg cm"^ s^S (64) 
which for Population I abundances {X = 0.7, Y = 0.28) gives 

f-. = 5.63 X lO^V^T'i^^ erg cm"^ s"\ (65) 

For the electron-electron bremsstrahlung, Svensson (1982) gives the following formula for 
Be < 1: 

fee = nlcrlmeC^af-^j^iU-37T^)el^^ 

X (1 + 1.161^ - 1.256*^/2) ergcm-^s'S (66) 

and for 9^ > 1- 

f- = nlcr'^meC^af24ee[ln{r]de) + 1.28] erg cm"^ s"\ (67) 

where = jm^c? is the classical radius of electron and r/ = exp(— 7^;) = 0.5616. Again, in the 
original formula there is a 5/4 instead of 1.28 (see NY). 



5.3. Synchrotron cooling 



Following NY we give formulae for the synchrotron emission of a relativistic Maxwellian 
distribution of electrons. These formulae are valid only for Q > \ but it is sufficient for applications 
to ADAFs (see NY). 

27r,^, ,di/3(r) 



where 



= 7^7^ O^XM, (69) 

2 ZTTrripC 
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and xm is the solution of the transcendental equation 



exp(i.88y9x^ j - 2.49 XiU ^ ^3^^2(1/^6) 



0.40 0.5316 \ 

7/6 17/12 I ' 



where radius r is in physical unit of cm. 



5.4. Compton cooling 



NY use the Dermer, Liang & Canfield (1991) prescription for the Compton energy 
enhancement factor ry: 



■q = 1 + 



I- PA 



1 



^ ,^ l-lnP/lnA' 

3^ 



= 1 + r?i + r?2 



(71) 



where 



P 
A 

m 
m 



hv 

= l-exp(-res), 

= i + We + ml, 

p{A-\) 
I- PA' 

= -l-lnP/lnA, 



(72) 



and Tes is the electron scattering optical depth. Note that 77, P, A here are different from the ones 
defined in previous sections. 

The Comptonization of the bremsstrahlung radiation is given by 



/br,c = ^mQbr 



3 SOeJ r?3 + l 



and the Comptonized synchrotron emission is 

/synch,C = /synch [^1 " ^ {Xc/&e)'"'] , 

where Xc = i^c/meC^. 



(73) 
(74) 
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5.5. Total radiative cooling 



In the optically thin case the total radiative cooling will be given by 



/ 2}{ -^br "I" /synch "I" /br,C "I" -^ynchjC 



whereas in the optically thick case one can use 

r-- 



3Ht 



(75) 



(76) 



where r = Tabs + '^es = 0-5{Kabs + «^es)S is the total optical depth. In the intermediate case one 
should solve the transfer equation to get reliable results. In the present context one can use the 
solution of the grey problem obtained by Hubeny (1990). NY give the following formula: 



r 



H 



37" 4o"^Z~"^ 

— + VS H ^— X (/j^j. + /gyn(,ij + /bi- c + /synch,c) 



where the relation: 



H 



''"abs 



■(/br + /synch + /br,C + fs- 



/^fj'J^4: br ' J synch Jhr,C J synch, C/ 



(77) 
(78) 



has been used. The opacity appearing in Eq. (^) is the Planck mean opacity. 
The radiation pressure can be expressed as 



f-H f 2 



(79) 



5.6. Thermal balance equations 

Since the viscous heating concerns the ions we shall write equation ( ^0|) in the form of two 
equations (NY): 

r = /-, (80) 

where /+ = F+/{2H), /^^^ = F^^/{2H), and f- are given by equations @) and (|^) 
respectively. 
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6. THE SONIC POINT AND THE BOUNDARY CONDITIONS 



Our system of equations, like its Newtonian and PW versions, possesses singular points. In 
this section we will use the vertically averaged quantities P and S. From equation (^) we have 

dlnS 1 dV 



where, 



dlnr V{l-V'^)dlnr 
^_ dln(AV2) ^ r{r-M) 



dlnr A 

The radial gradient of P in equation (36) can be derived from equation (^) and (46) as 



dlnP _ 2(4-3/3*)(ilnr 3/3*-ldlnS {1 - P*) dlnj'^g 
dlnr ~ 1 + dlnr ~^ 1 + /?* dlnr I + (3* dlnr 
Thus, by combining the radial motion equation (^) and the energy equation (^0|), one gets 

dlnV N 



dlnr D 



where 



D = Cci 



1 + 



2(1 -/?*) 



N 



-A - Bci 



2(1-/3* 



1 + /3* 



1 + 13* 

(Fg - l)47rr2 
(1 + I3*)M 



(81) 
(82) 

(83) 

(84) 
(85) 



(FH 



(1-/3*) 
(1 +/?*)' 



dl^g 2v(^) [dv(^)/dr^ 



A is defined by Eq. (|37|) and 



c^ = ^ 



1 + 2 



dlnr 



1 



1 



The equation for the temperature is then 



dlnr , ^ N , ^ 



1 + /3* 

27rr2(P+ - P^ 



(86) 



(87) 



McyT 



There are two singular points of these equations: horizon r = ru^ defined by A(r//) = 0, and 
sonic point r = rg defined by the condition 



V\rs) 



Cirs)clirs) 



l + 2c2(l-/3*)/(l + ^*) 



(89) 



(i) The horizon. At the horizon A = and therefore ;S = oo, ^ = oo. However, one may 
demonstrate that J\f = 1 and A/B = —1. Since on the horizon V'^{rj{) = 1, it follows that at the 
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horizon dV/dr is non-singular. For /?* 7^ 1 the gradients dTi/dr, dP/dr and dT/dr are singular 
only through the infinit boost required by the hydrostatic equilibrium equation (44). Note, that 
from equation ( p4[ ) it follows that at the horizon = 0. There are no extra regularity conditions 
at the horizon. 

(ii) The sonic point. For large radii the flow is subsonic, < c^, and as we have just seen, 
it crosses the horizon with the light speed, 1 = > c^. Thus, somewhere in the flow it must 
be = c^, i.e. the flow must pass through the sonic point. Derivatives dVjdr and dT/dr are 
non-singular at the sonic point only if, 

N = D = {). (90) 

This gives an extra algebraic condition which must be satisfied by an acceptable, regular at the 
sonic point, global solution of the slim disk equations. The regularity condition at the sonic point 
was the main source of numerical difficulties in solving the PW slim disk equations: only for one 
particular choice of the eigenvalue Cq the condition (pOD could be satisfied. 



One may demonstrate that in the limit of spherical accretion ($7 = 0, AA = 1) onto a 
non-rotating black hole (a = 0), the condition (]90| ) becomes (for (3* = 1), 



y\rs) = , ^ ^ (91) 
2[ts-\m) 



which in the notation of Shapiro h, Teukolski (1983) takes the form 

M 
2^ 



u^rs) = 1^, (92) 



that is identical with their formula (G.17). 

On the other hand using Eq. (|4^ ) in the form 

= clQ (93) 

and including into the Eq. (^) the dH/dr term: 

e.-l(lo,A'/^^).-^-V (.4) 

where M = dlogH/dlogr one obtains from Eqs. ( |8^ , ^) (for (3* = 1) the following condition for 
the velocity at the sonic point: 

V{rs) = (95) 
i + i 1 

which is the same condition as the one obtained by Narayan, Kato & Honma (1996). 

(iii) The boundary conditions. The set of equations describing the flow is of the fifth order with 
four independent radial first derivatives (say, dQ/dr, dV/dr, dTi/dr, and dT/dr) and one second 
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derivative {(PQ/dr'^). Five integration constants are required. The equation of mass conservation 
is integrated analytically, and this fixes one integration constant, M, which is a free parameter. 
Furthermore, since the angular momentum equation has been analytically integrated from the 
horizon, an (unknown) value of the specific angular momentum at the horizon, Cq, appears in the 
angular momentum equation. This constant however is not a free parameter and is determined 
such that the flow is transonic. The two regularity conditions {N = 0, D = 0) at the sonic point 
provide one condition for the sonic point location, and one extra condition for the five integration 
constants. Three boundary conditions are assumed at the outer boundary. There are no other 
boundary assumed conditions, in particular no other boundary conditions at the horizon. The 
solution is determined by M, a, and M, wheares Cq, is an eigenvalue that should be determined. 
In the numerical procedure, we use the self-similar solution (of Q, V, and T) as the outer boundary 
conditions, and by varying Co, we find its particular value (the eigenvalue of the problem) that 
fulfils the regularity conditions {N = 0, D = 0) (for detail see Chen & Taam 1993). 



7. NUMERICAL SOLUTIONS 

The set of Kerr black hole disk equations was solved by using the method described in Chen 
& Taam (1993). This method has been modified to take into account the terms introduced by 
general relativistic eff^ects. Here we focus only on the case of advection dominated accretion 
flows (ADAFs) and study the effects of the rotation of the black hole. Specifically, we assume 
the gas is optically thin (r = 0) and consider the /?* = 1 case.. The local radiative cooling 
only through non-relativistic bremsstrahlung emission (Eq. (|65|). The corresponding solution in 
pseudo-Newtonian potential approximation was obtained by Chen, Abramowicz, &: Lasota (1996) 
and Narayan, Kato, & Honma (1996). 

The mass accretion rate is measured in the Eddington limit which is deflned as 

Me = inGM/icKes), i^es = 0.34. (96) 

and the units for r and angular momentum (such as C, a) are tq and M respectively. 

We present in Figure 1 examples of solutions with three different rotation rates of the black 
hole, a/M = (solid- line) , 0.5 (dotted- line), and 0.99 (dashed- line) . The black hole mass is 
M = IOMq, the viscosity parameter is a = 0.1 and the accretion rate is M /Me = 10^^. The 
pseudo-Newtonian solution is also plotted for comparison (heavy dots). One can see that the 
pseudo-Newtonian solution is an excellent approximation of the solution for a Schwarzschild black 
hole (a = 0). As the black hole rotates faster, the sonic point moves inwards, mainly due to 
the increase of the sound speed. This is due to the deepening of the potential well depth with 
increasing black hole rotation since the ADAF temperature is close to the virial one. The specific 
angular momentum lost to the black hole is smaller for faster rotating black hole which is also the 
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effect of the hole becoming smaller with increasing rotation. At the black hole sm'face the pressure 
goes to infinity due to the infinite boost required by Eq. (^). 



8. CONCLUSION 

We have derived a complete set of equations describing accretion flows with non-zero angular 
momentum around a rotating black hole (non-Keplerian accretion disks). We used the Novikov 
and Thorne form of the Kerr metric, which is accurate up to the (H/r)^ terms. Since ADAFs can 
be heated up to the virial temperature and therefore H/r ^ 1, the slim disk approximation may 
appear to be not adequate. However, it is known from the study of ADAFs in the Newtonian case 
(Narayan & Yi 1995a) that the quasi-spherical structure of ADAFs can be well represented by 
suitably interpreted vertically averaged quantities. We don't expect general relativistic effects to 
modify this conclusion. 

We obtained several test solutions which show that ADAFs in a Kerr space time have the 
same basic properties as the pseudo-Newtonian ADAFs. The rotation of the hole allows however 
for higher temperatures and pressures which should modify the observed emission from ADAFs. 
This property as well as the influence of the general relativistic effects on light propagation will be 
studied in a future work. 

MAA acknowledges a partial support from Nordita through the Nordic Project on Non-Linear 
phenomena in accretion disks around black holes. 
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Fig. 1. — The radial structure of tlic Mach number (M), the pressure (p in cgs unit), the angular 
momentum (C in unit of M), and the sound speed (c^ in unit of the speed of light). Here the 
mass of the black hole is lOM©, a = 0.1 and M /Me = 10~^. The solid, dotted, and dashed lines 
represent the cases of a/M = 0, 0.5, and 0.99 respectively. The heavy dots represent solutions 
obtained with the pseudo-Newtonian potential. These solutions are excellent approximation to the 
solutions representing the Schwarzschild black hole flows (in the case a = 0). The corresponding 
Keplerian angular momenta of test particles around Kerr black holes are also shown for comparison 
(the thin lines). 



